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Magnetosphere-Ionosphere Coupling Through E'-region 
Turbulence: Anomalous Conductivities and Frictional Heating 

Y. S. Dimant and M. M. Oppenheim 

Abstract. 

Global magnetospheric MHD codes using ionospheric conductances based on laminar 
models systematically overestimate the cross-polar cap potential during storm time by 
up to a factor of two. At these times, strong DC electric fields penetrate to the E re- 
gion and drive plasma instabilities that create turbulence. This plasma density turbu- 
lence induces non-linear currents, while associated electrostatic field fluctuations result 
in strong anomalous electron heating. These two effects will increase the global ionospheric 
conductance. Based on the theory of non-linear currents developed in the companion pa- 
per, this paper derives the correction factors describing turbulent conductivities and cal- 
culates turbulent frictional heating rates. Estimates show that during strong geomag- 
netic storms the inclusion of anomalous conductivity can double the total Pedersen con- 
ductance. This may help explain the overestimation of the cross-polar cap potentials by 
existing MHD codes. The turbulent conductivities and frictional heating presented in this 
paper should be included in global magnetospheric codes developed for predictive mod- 
eling of space weather. 



1. Introduction 

In the companion paper [Dimant and Oppenheim, 
2011], we have developed a theoretical description of 
magnetosphere-ionosphere (MI) coupling through electro- 
static plasma turbulence in the lower ionosphere in the re- 
gion where field-aligned magnetospheric currents close and 
dissipate energy. In global MHD computer codes intended 
for predictive modeling of space weather, the entire iono- 
sphere plays the role of the inner boundary condition. The 
ionospheric conductances employed in these codes are usu- 
ally based on simple laminar models of ionospheric plasma 
modified by the effects of precipitating particles. These 
global magnetospheric MHD codes with normal conduc- 
tances often overestimate the cross-polar cap potential dur- 
ing magnetic (sub)storms up to a factor of two [e.g., Winglee 
et al, 1997; Raeder et al, 1998, 2001; Siscoe et al, 2002; 
Ober et al, 2003; Merkin et al, 2005a, b; Guild et al, 2008; 
Merkin et al., 2007; Wang et al, 2008]. At the same time, 
strong DC electric field penetrates to the E region, roughly 
between 90 and 130 km of altitude, and drives plasma in- 
stabilities. The .©-region instabilities create turbulence that 
consists of density perturbations coupled to electric field 
modulations. This turbulence causes anomalous conductiv- 
ity which could account for the discrepancy. 

Anomalous conductivity manifests itself in two ways. One 
is associated with average anomalous heating by turbulent 
electric fields, and the other is due to a net non-linear cur- 
rent (NC) formed by plasma density turbulence. While the 
former effect has been considered before in detail (see refer- 
ences below), the latter has only been discussed but never 
quantitatively investigated with application to ionospheric 
conductivity. This is the main objective of this paper. In 
the two following paragraphs, we outline each of these ef- 
fects. 

Small-scale fluctuations generated by the ©-region in- 
stabilities can cause enormous anomalous electron heating 
(AEH) [Schlegel and St. -Maurice, 1981; Providakes et al., 
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1988; Stauning and Olesen, 1989; St. -Maurice et al, 1990; 
Williams et al., 1992; Foster and Erickson, 2000; Bahci- 
van, 2007], largely because the turbulent electrostatic field, 
8E = — V(5<E>, has a small component parallel to the geo- 
magnetic field Bo [St. -Maurice and Laher, 1985; Providakes 
et al, 1988; Dimant and Milikh, 2003; Milikh and Dimant, 
2003; Bahcivan et al, 2006]. Anomalous electron heating 
directly affects the temperature-dependent electron-neutral 
collision frequency and, hence, the electron part of the Ped- 
ersen conductivity. This part, however, is usually small com- 
pared to the dominant electron Hall and ion Pedersen con- 
ductivities. At the same time, AEH causes a gradual ele- 
vation of the mean plasma density within the anomalously 
heated regions through the thermal reduction of the plasma 
recombination rate [Gurevich, 1978; St. -Maurice, 1990; Di- 
mant and Milikh, 2003; Milikh et al, 2006]. The AEH- 
induced plasma density elevations increase all conductivities 
in proportion by as much as a factor of two. This mecha- 
nism requires tens of seconds or even minutes because of 
the slow development of the ionization-recombination equi- 
librium. If Eo changes faster than the characteristic recom- 
bination timescale then its time-averaged effect on density 
will be smoothed over field variations. 

Also, low-frequency turbulence in the compressible iono- 
spheric plasma can directly modify local ionospheric conduc- 
tivities via a wave-induced non-linear current (NC) associ- 
ated with plasma density irregularities [Rogister and Jamin, 
1975; Oppenheim, 1997; Buchert et al, 2006]. The physical 
nature of NC has been explained in the companion paper 
[Dimant and Oppenheim, 2011]. While the rms turbulent 
field {SE 2 } 1 / 2 is comparable to Eo (the angular brackets de- 
note spatial and temporal averaging), the NC is proportional 
to the density perturbations that, in saturated FB turbu- 
lence, may reach tens percent at most. As a result, the total 
NC, considered as a plasma response to the external electric 
field Eo, amounts to only a fraction of the regular electrojet 
Hall current. However, in most of the electrojet the NC is 
directed largely parallel to Eo, so that it will increase signif- 
icantly the smaller Pedersen conductivity. This is critically 
important because it is the Pedersen conductivity that al- 
lows the field-aligned currents to close and dissipate energy. 
The combined effect of the NC and AEH makes the iono- 
sphere less resistive. These anomalous conductances may, at 
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least partially, account for systematic overestimates of the 
total cross-polar cap potential in global MHD models that 
employ laminar conductivities. 

In this paper, using expressions for the NC obtained in 
Dimant and Oppenheim [2011], we quantify a feedback of 
developed B-region turbulence on the global behavior of 
the magnetosphere by calculating the corresponding mod- 
ifications of local conductivities. In Sect. 2, we do this 
in general terms of given spectra of density irregularities. 
In Sect. 3, invoking the appendix, we employ a heuris- 
tic model of non-linearly saturated turbulence [Dimant and 
Milikh, 2003] which allows us to easily estimate the turbu- 
lent conductivities. The anomalous heating-induced effect 
manifests itself via the mean density variations so that it 
contributes proportionally into both the laminar and turbu- 
lent conductivities. In Sect. 4, we calculate the turbulent 
frictional heating sources for possible inclusion into global 
ionosphere-thermosphere models. In Sect. 5, we demon- 
strate that all anomalous effects combined may nearly dou- 
ble the undisturbed Pedersen conductance and hence they 
should be taken into account when applying inner boundary 
conditions in global MHD codes employed for space weather 
modeling. In the appendix, using a novel and compact for- 
malism, we derive the general fluid-model dispersion relation 
for waves in arbitrarily magnetized linearly unstable plasma. 

2. Turbulent Ionospheric Conductivity: 
Quasilinear Approximation 

In this section, we calculate the turbulent ionospheric 
conductivity associated with the B-region non-linear cur- 
rent, using the general quasilinear expressions for the NC 
obtained in the companion paper [Dimant and Oppenheim, 
2011]. This quasilinear approach is justified by the relatively 
small plasma density perturbations, even though electro- 
static field modulations can be comparable to the driving 
electric field [Dimant and Milikh, 2003; Oppenheim et al, 
2011]. 

The turbulent conductivity tensor can be constructed 
similar to the laminar conductivity tensor, 
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where the general parallel, Pedersen, and Hall conductivities 
in a Cartesian system £1,2,3 with the axis £3 directed along 
B are given by [e.g., Kelley, 2009] 
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We can find o using quasilinear Eq. (42) from Dimant 
and Oppenheim [2011] obtained for arbitrarily magnetized 
particles, 
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where, according to Eq. (24) from Dimant and Oppenheim 
[2011], 



foil = " 
Uo± = 



(« e + Hi) Eo\\ 
B 

{K e + Hi) [{Ke - 



{Eo X b) - (1 + KiKe) E ±] 



(l + K?)(l + /e?)B 



(6) 



3L\\ ■ E 0\\ _ {Ke + Ki 
jL± ■ E 0± 

El 



B 

{Ke + Ki){l + KiK e )ne 



l _ 3l± • {E x b) 

a H = S2 



{l + Kl){l + Kj)B' 



(2a) 
(2b) 
(2c) 



are the parallel and perpendicular to Bo components of the 
laminar relative velocity between the undisturbed electron 
and ion streams; U = VeO — Vio- In the quasilinear approach, 
the total conductivity, o , is determined by merely adding 
the laminar and turbulent conductivities, 
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respectively. Here n is the density of the quasi-neutral 
plasma; Jl±,\\ are the perpendicular and parallel to Bo com- 
ponents of the laminar current density jl\ k 3 = Qs/v s 
are the magnetization parameters for particles of the s-kind 
(s = e, i) with Q, s = eB jm s and v an being the corresponding 
gyro-frequency and mean collision frequency with neutrals; 
e is the (positive) elementary charge; m a are the particle-s 
masses. These conductivities neglect Coulomb collisions, as 
compared to electron and ion collisions with neutrals, so that 
they are largely applicable to altitudes below the .F-region 
ionosphere. 

Similarly to Eq. (1), the turbulent condu ctivi ty tensor 
determined by the non- linear current, J^ = ■ Eq, can 



where each of them is proportional to the same plasma den- 
sity no. Density increases caused by the AEH- induced re- 
duction of the recombination rate are automatically included 
if one allows for the corresponding temperature modifica- 
tions of no(T'e). 

Global MHD codes employed for magnetosphere mod- 
eling assume equipotential magnetic field lines, E u = 0, 
across the ionosphere due to high electron mobility along 
Bo- In the E/D regions, this approximation is valid above 
the 80 km altitude where k e 3> 1,k», i.e., almost in the entire 
electrojet. Hence, only conductivities perpendicular to Bo 
matter, while the parallel currents can be determined using 
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the continuity of the total current density, Vii = — V ±-j±. 
In this approach, the high parallel electron conductivity is 
not invoked unless one needs to estimate the tiny parallel 
field, S ||. 

Integration of the quas i-neutral current conservation 

equation, V-J tot = -V-(a tot -V<I>o) = 0, along nearly vertical 
equipotential magnetic field lines yields a 2-D second-order 

differential relation for the potential, Vi fs tot Vi$ol = 

where E tot = E L + S NC ^ is the total height-integrated iono- 
spheric conductance tensor and ju is the parallel current 
density on top of the conducting ionosphere. This height- 
integrated Ohm's law serves as an approximate inner bound- 
ary condition for the global MHD codes [e.g., Merkin et al., 
2005a]. 

With neglect of E \\, E ± ~ Eq = — V<&o, while assum- 
ing K e ^> 1 > Ki, the normal conductivity tensor given by 
Eq. (2) reduces to 
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Under the same conditions, while additionally presuming 
field-aligned irregularities, fcy <C k± « k, Eqs. (5) and (6) 
simplify to 
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where the last approximation of Eq. (10) represents the sim- 
plest interpolation between lower E/D-region altitudes with 
Ki < 1, K e > 1 > KiK e and higher .©-region altitudes with 
Ki ~ 1, KiK e 3> 1. Quasilinear Eq. (9) applies to arbitrary 
spectra of iJ-region turbulence, provided the density per- 
turbations are reasonably small (see the discussion below 
related to Fig. 1). In the next section, we employ a spe- 
cific model of non-linearly saturated turbulence to calculate 

NC 



3-D PIC simulations [Oppenheim et al., 2011]. A significant 
advantage of employing the HMT is that this model enables 
one to estimate the turbulent conductivities in terms of the 
rms density irregularities, without any knowledge of their 
detailed k, o;-spectrum. 

The HMT consists of two major heuristic assumptions 
of developed turbulence during the non-linearly saturated 
stage of the FB instability [Dimant and Milikh, 2003]. First, 
the model assumes the effective values of the major perpen- 
dicular component of the turbulent electric field, SEt ± -L B. 
Second, it assumes the effective aspect angles determined by 
effective k«/k±. Both these assumptions refer to a modified 
FB-instability threshold field within the developed turbu- 
lence, where the combined linear growth/damping rate 7g 
equals zero. This modified threshold field, hereinafter re- 
ferred to as merely the threshold field, is determined by 
the same expression as the actual threshold of instability 
excitation where the plasma temperatures are replaced by 
the elevated temperatures due to anomalous heating. We 
determine this threshold field using general expressions for 
the Farley-Buneman linear growth rate obtained in the ap- 
pendix. 

To calculate the threshold electric field necessary to main- 
tain turbulence, -E™h?> we can write 7g in the limit of 

Ke > 1 > Ki as 
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where xg i s t ne an gle between k and Uo [Dimant and Op- 
penheim, 2004], C's = (T e -Hi) 1 / 2 /mi is the isothermal ion- 
acoustic speed, Co = | Uo \, 
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following from Eq. (10), we can rewrite Eq. (11) as 
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where i?xhr is the minimum FB-instability threshold elec- 
tric field at a given altitude for optimally directed waves 
with k± || Uo and k\\ = (i.e., for = and ipg = tp±), 



3. Turbulent Ionospheric Conductivity: 
Heuristic Model of Turbulence 

The current state of ©-region instability theory does not 
give us accurate spectra of density irregularities 8n^ ^ as 
functions of the external electric field and ionospheric pa- 
rameters. In order to estimate the turbulent conductivities 
we are forced to use simplified models of non-linearly satu- 
rated turbulence. Here we employ a heuristic model of tur- 
bulence (HMT) developed previously for qualitative expla- 
nation of AEH by Dimant and Milikh [2003] . The HMT pro- 
vides approximate rms values of the coupled turbulent elec- 
tric field and density irregularities. This analytical model 
based on simple physical reasoning has been successfully 
tested by quantitative comparisons with AEH observations 
[Milikh and Dimant, 2003; Milikh et al, 2006] and our recent 
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Here 
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(16) 

is the absolute threshold-field minimum which can only be 
reached at optimum altitudes for the FB instability excita- 
tion where the conditions of k- C 1 and ip± — Oq/k 2 <§C 1 
overlap (at high latitudes, this occurs at 100-105 km al- 
titude [Dimant and Oppenheim, 2004, Fig. 2]; O 2 , = 
m e u e /(miVi) ~ 1.8 x 10~ 4 ). Be advised that notations in 
this paper slightly differ from those in Dimant and Milikh 
[2003, Eqs. (14), (15) and others]. Namely, Eo in this paper 
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corresponds to Ec in Dimant and Milikh [2003], while our 
^Thr corresponds to Eq and £?xhr corresponds to Bxhr in 
that paper. 

The first heuristic assumption specifies the rms perpen- 
dicular turbulent electric field, 



(SEl) 



ai (E -E?£) 2 , 
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where ai is a dimensionless factor of order unity [Dimant 
and Milikh, 2003, Eq. (25)]. The logic behind Eq. (17) 
is that in the non-linearly saturated state the major elec- 
tron non-linearity oc (8E± x Sn) balances, on average, the 
linear instability growth for optimally directed waves. For 
-Bo > E^, the rms turbulent field, {SE\) l/2 , is similar in 
magnitude to Eq, while near the minimum FB instability 
threshold, E « E%£, it reduces linearly with (E - E%£). 
This heuristic assumption makes no distinction between the 
2-D and 3-D cases, implying that the rms perpendicular tur- 
bulent fields in both cases are approximately the same. This 
has been confirmed by our recent fully kinetic simulations 
[Oppenheim et al, 2011]. 

The other heuristic assumption distinguishes the 3-D case 
from the purely 2-D one by quantifying the effective magni- 
tude of fcii . In B-region turbulence, fey is always much less 
than k± « k, but even small k« matter because for highly 
mobile electrons with K e ^> 1 they significantly modify the 
parameter ipg defined by Eq. (12a) and, hence, the linear 
instability threshold. Also, it is the parallel turbulent field, 
5E\\ oc fcii , that largely cause AEH. The idea of the second 
HMT assumption is that the effective values of k\\ in devel- 
oped turbulence, on average, settle the system on the margin 
of linear stability where 7g = 0. Given the optimum value 
of the flow angle for the pure FB instability, the marginal 
linear stability yields 
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(cf. Dimant and Milikh [2003, Eq. (26)] with our ipf corre- 
sponding to (1 + k 2 )Vwx in Dimant and Milikh [2003] (the 
conditions in the two papers look formally different, but they 
are essentially the same because of the negligible difference 
~ ip±K- = 6q ~ 1.8 x 10 -4 ). Equations (27)-(29) from Di- 
mant and Milikh [2003] yield specific values for the effective 
rms values of the parallel turbulent field in terms of a second 
unknown constant of order unity, «2 • That was important 
for determining the average heating source responsible for 
AEH but is not required for our current purposes. Equa- 
tions (17) and (18) is all we need from HMT to quantify 
the rms density perturbations and hence the corresponding 
non-linear current. 

In the quasilinear approximation, the turbulent electric 
field spectrum, SE? — —ikS&g^, is proportional to the 
density irregularity spectrum Sn^ ^ . In our current limit of 
K e > 1 > Ki, Eq. (28) from Dimant and Oppenheim [2011] 
yields 
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where the proportionality coefficient between SE^ and 
Sng u depends on the wavevector direction via cos \% ~ 



k-U~o/(kUo) and ijjg, but is independent of its absolute value, 
k. This fact allows us to directly relate the rms value of 
the entire turbulent field, {SE 2 ± ) 1/2 = (E*Wo l 5 %J 2 ) 1/2 , 
to the rms of the total density fluctuations, {Sn 2 ) 1 / 2 = 
(Ej \& n k J 2 ) 1 ^ 2 - Indeed, if we assume in accord with 
simulations [e.g., Oppenheim and Dimant, 2004; Oppenheim 



et al., 2008, 2011] that most of the fex-spectrum is concen- 
trated within a narrow angular sector around a preferred 
wavevector direction 
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where the flow angle \^ with the positive sign of sinx? 1 



corresponds to the tilt in the (— £ , )-direction [Dimant and 
Oppenheim, 2004], then, setting for the entire spectrum 
fc/fe ~ (k/k) pr , tpj: ~ -i/>™ and using Eq. (13), we obtain 
from Eq. (19) an approximate relation 
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According to Eq. (15), we have 
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where B^hr is defined by Eq. (16). 

Figure 1 based on Eqs. (21) and (22) with ai = 1 and 
Xf = [Dimant and Milikh, 2003] shows that if E /E^ T 
is large enough then relative rms density perturbations may 
even exceed unity. Such non-physical occasions break the 
validity of our quasilinear approximation. We should bear 
in mind, however, that large values of Eq automatically 
lead to strong electron and ion temperature elevations due 
to combined regular and anomalous heating. This raises 
E tL « ( r <= + Tt) 1/2 and prevents Eq/E^ from reach- 
ing too large values (we discuss this issue in more detail 
in the next section). For example, in the extreme case of 
i?o = 160mV/m, the electron temperature reaches above 
4000 K [Bahcivan, 2007] at an altitude around 110 km corre- 
sponding to m ~ 0.3 [Dimant and Oppenheim, 2004, Fig. 2], 
not counting the simultaneously increasing ion temperature. 
According to Eq. (16), this raises B™ r above 50mV/m, 




Figure 1. Relative density perturbations, {Sn 2 ) 1 / 2 /no, 
vs. the ion magnetization parameter m = Qi/vi for dif- 
ferent values of Eo/E™^ (shown near the curves). 
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making Eq/E^ t < 3. We believe that this extreme value 
of -Eo/^xhr imposes the top restriction on possible density 
perturbations, keeping them below 50%. Less extreme but 
more typical density perturbations are expected to be largely 
within the reliable limits of the quasilinear approach. 

Applying k/k ~ (k/k) pT , while replacing (1 + ^j?) accord- 
ing to Eq. (18) and (Sn 2 ) 1 ^ 2 according to Eq. (21), we use 
Eqs. (9) and (10) to express the turbulent conductivities in 
terms of the corresponding laminar conductivities given by 
Eq. (8) as 
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These relations are only applicable for Eo > -E™hr an d K i < 
1, otherwise apg = 0. Furthermore, according to Eqs. (15) 

and (16), the threshold field E^£ °c ^xhr « ( T <= + Ti) 1/2 de- 
pends strongly on the electron and ion temperatures. This 
means that the self-consistent inclusion of the turbulent con- 
ductivity in the total conductivity tensor requires simultane- 
ously including AEH [Dimant and Milikh, 2003; Milikh and 
Dimant, 2003], otherwise the effect of NC-induced turbulent 
conductivity will be exaggerated dramatically. 

The AEH-induced increase in the plasma density due to 
the temperature reduction of the recombination rate [Di- 
mant and Milikh, 2003; Milikh et ai, 2006] will raise all val- 
ues in the conductivity tensor in proportion to the increased 
plasma density. This slowly-developing effect, however, can 
slightly decrease if the strong external electric field Eo varies 
too fast [Codrescu et al, 1995, 2000; Matsuo and Richmond, 
2008; Cosgrove et al, 2009; Cosgrove and Codrescu, 2009]. 



4. Turbulent Prictional Heating in Global 
Modeling 

An accurate description of strong electric-field pertur- 
bations, such as storms, substorms, and sub-auroral polar 
streams, requires including macroscopic effects of i5-region 
turbulence into global computer models, like the Coupled 
Magnetosphere Ionosphere Thermosphere (CMIT) model 
[Wiltberger et ai, 2004; Wang et al, 2004]. This model 
has been created by combining the magnetosphere MHD 
solver LFM [Lyon et al, 2004] and ionosphere-thermosphere 
NCAR solver TIEGCM [Roble and Ridley, 1994; Wang 
et al, 1999]. The ring-current model RCM [Toffoletto et al, 
2003] has been added recently, as well as TIEGCM has been 
extended to cover the mesosphere (TIMEGCM). Currently, 
the CMIT-RCM model includes ionospheric processes asso- 
ciated only with laminar conductivities. The turbulence- 
induced non-linear current (NC) can be incorporated by 
adding turbulent conductivities given by Eq (23). Self- 
consistency requires including also anomalous plasma heat- 
ing caused by turbulent fields. The AEH increases (through 
the reduced recombination rate) the plasma density and, 
hence, all conductivities in proportion. However, the anoma- 
lous temperature elevations exert a negative feedback on the 
saturated level of the plasma turbulence, reducing the NC 
and its non-linear Pedersen conductivity. 

Effects of i?-region turbulence can be included directly in 
TIMEGCM and other codes that that model high-latitude 



neutral and plasma dynamics and chemical reactions. Given 
accurate heating sources, TIMEGCM includes all ioniza- 
tion/recombination and collisional cooling processes auto- 
matically. To account for temperature inputs caused by 
both laminar and turbulent fields, the energy balance equa- 
tion for each species should include the corresponding source 
of frictional heating. Currently, TIMEGCM includes the 
laminar ion and neutral frictional heating but neglects the 
electron one. This is a reasonable approximation for al- 
titudes above 100 km under quiet ionospheric conditions. 
However, during strong electric-field events, the i5-region 
turbulence via AEH dramatically raises the electron tem- 
perature, meaning that the total electron frictional heating 
should also be included. For the convection field well above 
the instability threshold field given by Eqs. (23), ion turbu- 
lent heating can also be appreciable (for Eo 3> E™^:, it is 
comparable to the ion laminar heating) . Accuracy of the en- 
tire global model requires including the corresponding neu- 
tral frictional heating as well. Note that highly anisotropic 
S-region turbulence also causes average momentum changes 
in the plasma. Plasma-neutral collisions in turn can transfer 
these changes to neutral particles. For the weakly ionized 
E region, however, relative momentum changes are much 
less important than relative energy changes, so that we will 
ignore the former and focus entirely on the latter. 

In the multi-component i?-region ionosphere where 
electron- neutral (e-n) and ion-neutral (i-n) collisions domi- 
nate, we can calculate the total laminar and turbulent fric- 
tional heating of electrons, H l e ot , ions, Hj ot , and neutrals, 
Hl? % , using a quasilinear approximation to quadratic accu- 
racy. In terms of the average turbulent perturbations, these 
are given by 

e v e [n (V 2 + (SV 2 )) + 2(V e0 ■ (SndV e ))] , 
^ rm + m n 

n 

x [n l0 {V 2 o + (SV 2 )) + 2{V i0 ■ {8n t 8V))\ , (24) 



rrijUin 
mi + rn n 



x [n M {V 2 + (SV 2 )) + 2(V ■ (SmSfy)] , 

with the summations taken over all possible collisions be- 
tween ions and neutrals, where indices i and n refer to sep- 
arate groups of ions or neutrals; mo = no. The physical 
meaning of various terms in the square brackets is explained 
in Dimant and Oppenheim [2011]. The expressions for H^°* 
explicitly take into account the fact that only m„/(m,+m„) 
of the partial i-n electric field heating, miVi„V 2 , goes to ions 
(//* ot ), while the remaining fraction goes directly to the col- 
liding neutrals (_ff^ ot ). The total electron collision frequency 
u e includes all electron-neutral (e-n) collisions, v e — v en ; 
but because m e <C m n , the contributions to from the 
e-n collisions are negligible. Equation (24) is written in the 
neutral frame of reference, presuming that all thermosphere 
components move with a common neutral-wind velocity. 

In a multi-fluid plasma, the laminar and turbulent veloci- 
ties, Vso and SV S , can be expressed in terms of the convection 
field Eo and spectral harmonics of the wave potential (Jc&g 
similarly to the two-fluid plasma [Dimant and Oppenheim, 
2011]. However, the first-order linear relations between the 
harmonics of the ion spectral densities, Sn^ w , and poten- 
tial, <5<3>g u , which are crucial for obtaining equations like 
Eq. (19), (43), can be rather complicated or cannot even 
be expressed in any closed analytical form. The reason is 
that the algebraic order of the underlying linear dispersion 
relation for the first-order wave frequency, cjg, in the gen- 
eral case equals the number of separate ion species. Even for 
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the two dominant iJ-rcgion ion species, NO+ and Oj, when 
the dispersion relation reduces to a quadratic equation, the 
first-order relations become cumbersome. These relations, 
however, simplify dramatically if all ion species have a com- 
mon magnetization parameter, m — Sli/i/i (i>i = z/; n ), so 
that all ion species respond to the fields equally. As a result, 
the multi- fluid relation between <5<3>r and Snr and vari- 

ous terms in n a0 (V 2 + (SV 2 )) + 2(U s o • (Sn s SV s )) (s = e, i) 
reduce to two-fluid Eqs. (49) and (51) from Dimant and 
Oppenheim [2011]. For K e j$> 1 > m, these become 

2 v e n El 2 ViUioKjEl 

Ve.n V e0 « — — — , nn m V i0 = j— — 2— ^, (25) 



B 2 



(1 + K 2 )B 2 ' 



v e n (SVe) 



miVin (l + nf) tp^jk ■ Up) 2 

" : ^ (i + ^) 2 fci 



<5n E 2 



n 



(26a 



^n i0 (<5Vf ) « ^n i0 (l + k?) ]T (l+'^°)2fe3_ 



n 



(26b) 



2en (l + «?) 



/Tie 



E 



(fc-g )(fc- t/ ) 



n 



(27a) 



2z/i(y i0 • {Smsvi)) 



2en i0 



E 



(k-E )(k-U ) 



2 
-L 



no 



(27b) 



where 5ng = <5 n i£ u is the electron density, V'g and Up 
are given by Eqs. (12) and (13), and (1 + k 2 )(1 + K 2 fc|/fci) 
emerging in the derivation of Eq. (26a) was approximated 
by ipg/ip± — K e KiVfe [Dimant and Oppenheim, 2011, see 
Eq. (34b) and text below it] . 

To calculate heating rates in terms of i?-region parame- 
ters, we need estimates of the turbulent density perturba- 
tions. Using the heuristic model of saturated turbulence, 
along with the approach that lead us from Eq. (9) via 
Eqs. (17), (18), and (21) to Eq. (23), we obtain 



(1 



k,u>^P 



5n k,c 



n 



Eq 

77mm 
^Thr 



(i + V±) - 1 



v (fc- t/p) 2 



n 



a^UEo-E^f cos 4 X f 

(l + K 2 ) 3 ^ 2 



ai KUEo~E^f cos 4 



(1 + K 2 ) 3 B 2 



E 



k,ui^0 



(k-E )(k-U ) Sr%u 
(1 + 4> u )k 2 ± n„ 



Q 1 K 2 (g^-£o) 2 cos^ X | r sin X g r 
(1 + k 2 ) 5 / 2 B 



(28) 



Equations (24)-(28) provide frictional heating terms for the 
potential inclusion into the corresponding energy-balance 
equations of TIMEGCM or similar ionosphere-thermosphere 
models. Equations (23)-(24) provide a reasonable estimate 
of turbulent conductivity and heating and should greatly 
improve modeling of the geospace environment during dis- 
turbed conditions. 



5. Discussion 

Figures 2-5 illustrate anomalous heating and give exam- 
ples of the anomalously modified Pedersen conductivity for 
two cases: an extreme field of 150 mV/m and a still large but 
more modest case of 80 mV/m. These figures are based on 
typical high-latitude ionospheric parameters and the same 
values of the major HMT parameter, ai = 1, as in Milikh 
and Dimant [2003]. Figure 2 shows the effect of AEH and 
normal ion heating. The combined temperature, T e + Tj, 
shows kinks at altitudes where AEH sharply disappears, 
while ion heating continues its steep rising. These altitudes 
are slightly below the ion magnetization boundary (m = 1) 
located at hi ~ 122 km. These kinks translate to those in 
the threshold electric field, E™£ oc (T e + T;) 1/2 , as seen in 
Fig. 3. Figure 3 (a) shows that above 100 km of altitude the 
threshold field increases dramatically due to plasma heating, 
keeping Ep/E™^ at a modest level of 2-3. The much weaker 
field of Eo = 80 mV/m (b) gives a slightly lower ratio be- 
cause the plasma heating in this case is also much lower. 
According to Fig. 1, such moderate ratio of Ep/E™^., even 
for the extreme case of Eo = 150 mV/m, keeps the total 
density fluctuation level below 50%, i.e., roughly within the 
framework of the quasilinear approach (see above). This 
also reduces the entire effect of NC on the conductivity, as 
seen from Fig. 4. 

Figure 4 shows a typical altitudinal profile of the undis- 
turbed conductivity (curve 0) along with two curves includ- 
ing the anomalous conductivity. The NC effect alone (curve 

1) nearly doubles the conductivity in a broad altitude range 
between 90 km and the ion magnetization boundary located 
at an altitude hi ~ 122 km, fairly close to the maximum 
of the normal Pedersen conductivity. The AEH-induced 
ionization-recombination effect doubles the total conductiv- 
ity once again (curve 2), but this occurs in a slightly nar- 
rower range restricted from above by the AEH upper bound- 
ary, ho ^ 119.5 km. Presuming that the lower altitudes lo- 
cated below the undisturbed conductivity maximum near hi 
include about a half of the entire height-integrated conduc- 
tance, we see that the NC-induced anomalous conductivity 
alone (curve 1) can contribute up to a half of the undisturbed 
conductance, while both anomalous effects combined (curve 

2) can nearly double the entire conductance. We should bear 
in mind, however, the caveats associated with slow evolution 
of ionization-recombination processes superposed on rapidly 
varying electric field (see the discussion above). 

Largely the same situation takes place for the smaller 
field of Eo = 80 mV/m (Fig. 5). Since in this case the 
particle-temperature elevation is more modest than that for 
Eo = 150 mV/m, the instability threshold field shown in 
Figure 3 (b) remains noticeably lower than that in Fig. 3 (a). 
As a result, the ratio Eo/E™^ changes much less dramati- 
cally compared to Bo = 150 mV/m. That is why in Fig. 5 
the NC-induced effect for Eo = 80 mV/m turns out to be 
nearly as strong as that in Fig. 4 for E = 150 mV/m, while 
the more fragile AEH-induced recombination effect is much 
less pronounced because of lower electron heating. 

Above we presumed a rather conservative value for the 
major HMT parameter, ai = 1, the same as in Dimant 
and Milikh [2003] ; Milikh and Dimant [2003] ; Merkin et al. 
[2005b]; Milikh et al. [2006]. Recent 2D and 3D PIC sim- 
ulations [Oppenheim et al, 2011] suggest that ai can ac- 
tually be closer to 1.5. This would make <Tp C /<7p roughly 
50% larger, although due to larger (Sn 2 ) 1 ^ 2 /no one may have 
stronger restrictions on the quasilinear approach. We should 
also note that our current heuristic model of developed tur- 
bulence does not include the ion thermal driving mechanism 
(ITM) [Dimant and Oppenheim, 2004; Oppenheim and Di- 
mant, 2004]. 
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Figure 2. Altitude dependence of the electron and 
ion temperatures, T e and T,, for two different values 
of Eo (shown in figures). The electron temperature in- 
cludes self-consistent AEH according to Milikh and Di- 
mant [2003]. The kinks in T e and T e + Ti near the 
119.5 km altitude (a) and 117.5 km (b) are caused by 
sharp disappearance of instability and, hence, AEH. 




Figure 3. Altitude dependence of the FB-instability 
threshold electric field, E^tr, Eqs. (15) and (16), ele- 
vated due to particle heating for the same values Eo as 
in Fig. 2. The kinks in oc (T e + T») 1/2 correspond 

to those in Fig. 2. 
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Figure 4. Altitude dependence of the total - nor- 
mal and anomalous - Pedersen conductivity for Eo = 
150 mV/m (in relative units). Curve shows the undis- 
turbed conductivity. Curve 1 includes additionally the 
NC-induced anomalous conductivity calculated accord- 
ing to Eq. (23a) with on = and xg = 0. Curve 2 shows 
the total Pedersen conductivity with the AEH-affected 
plasma density elevated according to the steady-state ion- 
ization recombination model by Milikh et al. [2006]. The 
NC-induced anomalous conductivity disappears above 
the ion magnetization boundary ~ 122 km; the AEH- 
recombination effect vanishes above the top boundary of 
anomalous heating ~ 119.5 km (see the corresponding 
kinks in Figs. 2 and 3). 
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Figure 5. The same as Fig. 4, but for E — 80 mV/m. 



That is why all our anomalous effects occur strictly below 
the ion magnetization boundary, Ki = 1. Possible inclusion 
of the ITM would play a two-fold role. On the one hand, this 
would expand the altitudinal range of anomalous conductiv- 
ity to at least a few kilometers higher [Dimant and Oppen- 
heim, 2004] and increase the total ionospheric conductance 
accordingly. A modest reduction in Ki due to anomalously 



heated ions might also help. On the other hand, near the ion 
magnetization boundary the preferred values of the modified 
flow angle x% deviate from xg = 0, also due to the ITM as 
explained in Dimant and Oppenheim [2004]. This deviation 
of x% would slightly reduce the xgi Ki-dependent factor in 
the RHS of Eq. (23a) and hence lead to a smaller Peder- 
sen conductivity compared to that for xj — 0. An accurate 
inclusion of the ITM to the heuristic model of saturated tur- 
bulence, however, would complicate the entire model. We 
plan to improve the treatment of the top electrojet alti- 
tudes in future by using both analytical theory and PIC 
simulations. At this moment, a better accuracy is of less 
importance than the mere fact that the instability-induced 
conductivity occupies roughly the entire lower half of the 
Pedersen conductive layer and can nearly double the whole 
conductance. Note also that polarized sporadic-i? clouds 
[Dimant and Oppenheim, 2010], or even ubiquitous meteor 
trails [Dtmant et al., 2009], can also make additional contri- 
butions to anomalous conductances. All these effects com- 
bined can, at least partially, explain why global MHD codes 
developed for predictive modeling of space weather which 
use normal conductances often overestimate the cross-polar 
cap potentials by close to a factor of two. 

6. Summary and Conclusions 

Plasma turbulence generated by S-region instabilities sig- 
nificantly modifies the ionospheric Peder sen c onductance by 

exciting a net non-linear current, J NC — ap^ ■ Eo- The mag- 
nitude of j^ can only be a fraction of the undisturbed Hall 
current, but can be comparable in magnitude to the Ped- 
ersen current and pointing in that direction. This creates 
an additional turbulent conductivity. Also, anomalous elec- 
tron heating, via reduced plasma recombination, increases 
the mean plasma density no and all, laminar and turbulent, 
conductances in proportion with no. Estimates based on a 
quasilinear theory and heuristic model of non-linearly sat- 
urated plasma turbulence show that the anomalous effects 
combined can nearly double the Pedersen conductance, as 
predicted by Eqs. (9) and (23) and illustrated by Figs. 4 
and 5. Such ionospheric response to the magnetospheric 
field Eo may efficiently reduce the high-latitude ionospheric 
resistance and decrease the cross-polar cap potential. This 
effect might explain, at least partially, why routine MHD 
simulations with the normal ionospheric conductances sys- 
tematically overestimate this potential by a factor of two. 
The anomalous effects on the conductances, as well as tur- 
bulent frictional heating, should be included in global MHD 
codes developed for space weather prediction. 

Appendix: General FB/GD Dispersion 
Relation for Arbitrarily Magnetized Plasmas 

This appendix develops the two-fluid linear theory of the 
Farley-Buneman (FB) and gradient drift (CD) instabilities 
for arbitrary magnetization parameters, n eti = £l e ^/v eti , 
that vary from small values at lower altitudes to large values 
at higher altitudes. This analysis has never been published 
but is relevant to the lower ionosphere and enables one to ac- 
curately estimate conventionally neglected terms that may 
matter in some cases. 

The approach described below allows one to obtain the 
fluid-model dispersion relation in a reasonably compact and 
physically clear way. Compared to Sect. 3.1 of Dtmant and 
Oppenheim [2011], here we additionally include the particle 
pressure, inertia, regular gradients in the background plasma 
density, and recombination. These factors are responsible 
for the FB and CD instability drivers, as well as for the sta- 
bilizing diffusion and ionization balance. For simplicity, we 
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do not include non-isothermal processes responsible for driv- 
ing of thermal instabilities [Dimant and Sudan, 1997; Kagan 
and Kelley, 2000; Dimant and Oppenheim, 2004] and ignore 
the stabilizing effect of a weak non-quasineutrality [Rosen- 
berg and Chow, 1998; Kovalev et al, 2008]. These factors 
can be added using essentially the same approach. 

The quasineutral two-fluid continuity equations including 
ionization and recombination balance are 



dm + V • (nV e ) = <9 t n + V • (nV t ) = Q- an, 



(29) 



where Q is the total source of ionization and a is the re- 
combination constant. Equation (29) includes the standard 
quasineutral relation V • (nU) = 0, where U = V e — Vi. 

The continuity equations for a stationary (dtno = 0) 
but inhomogeneous background density, no, with the cor- 
responding undisturbed fluid velocities, 14, ;o, are 



V • (noKo) = V • (n Uio) 



2 

aon . 



(30) 



From the left equality, the divergence of the relative velocity, 
Uo = V e o — Vio, can be expressed in terms of the background 
density gradient as 



Now we express 8V e ; g in terms of <5ng and <5$g w , 8E^ = 

—ikS&j: w , using the momentum equations for the individual 
electron and ion mean fluid velocities in the neutral frame, 



d e U £ 



-e(E + V e x B) - V (nTe) - m e v e V e , (38a) 
dr n 



at n 



— nnViVi, 



(38b) 



where d e< i/dt = dt + V e ,i ■ V are the full derivatives for 
both electrons and ions. For low-frequency E/D-region 
plasma processes, the electron inertia described by the LHS 
of Eq. (38a), unlike the ion inertia in Eq. (38b), is usually 
neglected, but we will keep it for completeness and symme- 
try. Neglecting in Eq. (38) temperature perturbations, we 
express the velocity perturbations as 



Gi 



k.cj 



e n 
Ti Sn^ 
e n 



(39a) 
(39b) 



V • Uo = — Uo ■ Vno/no. 



(31) where the vector-functions G e 



In what follows, we presume that both Uo and no are known 
functions of spatial coordinates that satisfy Eq. (31). Then, 
to express the zero-order drift velocities of electrons and ions 
in terms of Uo, we can use Eqs. (26) and (27) from Dimant 
and Oppenheim [2011] to obtain 



oil 



KiU \\ 



(32) 



G e 



m e fc|| 
B~ 



ik e (k± — k e k x b) 

iki(k± + kik x b) 
~ (1 + k?)B ' 



(40a) 
(40b) 



describe anisotropic two-fluid responses to the harmonic per- 
turbations of the potential and pressures combined. Here, 
in accord with Eq. (45b), 



Ke[Up± - m(Ua X b)} 



v i0± = - 



Ki[U ± + K e (Up X b)] 



(33a) 
(33b) 



Note that these expressions are only valid when neglecting 
for the background plasma pressure gradients and gravity, 
comparably small effects for almost all ionospheric condi- 
tions. 

We will now consider Fourier harmonics of wave per- 
turbations, Snt, 8E^, 5V 4 g oc exp[i(k • r — wrt)], where 
the wavevector k and the complex linear wave frequency, 
cjg = u)£ + i7g, are locally defined. Here cjg is the real wave 
frequency, while 7g = —i Im o;g is the wave total growth or 
damping rate. 

Linearizing Eq. (30) with respect to 8n^, SV e ^ and in- 
troducing two shifted complex wave frequencies, 



)(«) 
k K 

nf =oj u -K- V i0 + 2ian = nf + q, 



(34a) 
(34b) 



with 

we obtain 



where 



K = k-iV, q = K-U , (35) 
= P-5V i% , (36) 
(37) 



k no k no 



Kno t> iVno 

p = — k . 

no n 



n e 



k 



(e) 



1 + 



iQ^ 



n t 



in 



1 + 



(41a) 
(41b) 



are modified magnetization ratios that include small con- 



tributions from the particle inertia; oL e ' 1 ' = Re f2- e ' lj — 

kik 

ojr — k ■ V e ,io- The small ion inertia contribution described 
by ifil'^/z/j is crucial for excitation of the FB instability. 
Now we substitute the expressions for SV e f g from 



Eq. (39) to Eq. (36), expressing nf = nf - q via 

q — k ■ Uo — zV • Uo = Uo ■ (p)* , where we used Eq. (31) and 
(35). This yields two independent linear relations between 

5n k and S ^k,c- 



4 e) + -(p-G e ) 

K p 



nf +q - 



(P- G t ) 



Snz , -* , 
Sn k 



no 



From these relations and Eq. (34b), we obtain two symmet- 
ric expressions for the coupled shifted frequencies, 



n^ = - 

k 



(?-G e )[g-ff-Gi)(r e + rQ/e] 

P ■ (Ge - Gi) 



(42a) 



n (0 = _ (p-GQ[g-(p-G e )(r e +r t )/ e ] 
% p-(Ge-Gi) 
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and the general relation between harmonics w and 8n^, Using Eqs. (12), we obtain 
q-p- {G e T e + GiTi)/e 5n e 



p ■ (G e - d) 



n 



(43) 



Using the definitions of Qr^' or fil 1 ' given by Eq. (34), we 



obtain the complex wave frequency in the neutral frame, 
Lo t = ni e) + K ■ V e0 - 2ian = + K ■ V i0 - 2ian , 

K k k 

{p-G e ){k-V i0 )-{p-Gi){k-V e0 ) 

(44) 



°k =UJ k+ m 



p-(G e ~G l ) 



p- (G e - Gi)e 

Equations (42) to (44) give the general expressions for the 
complex wave frequencies. To find the combined linear 
growth of instabilities or wave dissipation, 7^ = —i Imwg, 
they should be further split into the real and imaginary 
parts. 

The approximate two-fluid description of local quasi- 
harmonic waves is valid provided the characteristic wavevec- 
tors and frequencies satisfy 



r- 1 ,|7 E l«Kl«^e, i . 



(45a) 
(45b) 



Here l e ,i are the typical mean free paths of the correspond- 
ing particles with respect to ion-neutral and electron-neutral 
collisions, p e ,i are the gyroradii (if the corresponding parti- 
cles are magnetized), and T and Lii± are typical temporal 
and spatial (parallel and perpendicular to Bo) scales of iono- 
spheric density variation. Equation (45a) should be satisfied 
separately for the parallel and perpendicular directions. As- 
suming these conditions, we will treat the wave pressure 
gradients, local gradients in the background density, and 
particle inertia as second-order effects with respect to the 
small parameters defined by Eq. (45). In this treatment, 
the first-order factors discussed in Dimant and Oppenheim 
[2011] define the real wave frequency, cjg, while all second- 
order factors combined define the total linear damping or 
growth rate, 7g. 

Under conditions specified by Eq. (45), the dominant real 
part of the wave frequency, w?, is determined by the first- 
order accuracy, u'^ « wgj , when neglecting the pressure gra- 
dients, particle inertia, regular gradients of the background 
plasma density, and recombination. To this accuracy, we 
have k « p « k, q « k ■ Uq, R e ,i = t e ,j, and hence 



G e || « G e l II = 



Gi|| w Gji|| 



iK K ku 



B ' 

lKik\\ 

B 



G c 



G e ix = 



iK e (k± — n e k x b) 
(46a) 

iKi(k± + Kik x b) 
(46b) 



p- G e « k ■ Gel 

p-Gi^k-Ga 



irKip k I 



(l + Ki)B 



i + + 4 



1 + (1 + K?) fcf 



(47) 



p-(G e -Gi) « k- (Gei - Ga) - {1 + K 2 ){1 + ^ B 

(48) 

After neglecting small terms proportional to (T e + T») and 
ano, Eq. (44) yields 



we ~ we, 



(fc-G e l)(fc-Vj )-(fc-Gil)(fc- Kq) 
fc- (Gel -Gil) 



and reduces to Eq. (37) from Dimant and Oppenheim [2011] 



with the corresponding real shifted frequencies, 
Refill = w c , - fe • Vt.io, 



(e,i) 



(fc- U )(k-G e 



kl k- (Gel -Gil) 



n<?> 
fci 



(fc-!7o)(fc-Gii) 

fc- (Gel -Gil) 



(49) 

given explicitly by the companion paper Eqs. (36) and (39) . 

Now we will develop the dispersion relationships to the 
second-order accuracy by taking into account all previously 
neglected factors: the wave pressure gradients, particle in- 
ertia, gradients of the zero-order plasma parameters, and 
recombination. All second-order terms in Eq. (44), repre- 
senting linear corrections to the first-order real frequency 
wr-p are purely imaginary, so that their linear combination 
determines the total linear growth/damping rate, jg. 

We start by discussing the two last terms in the RHS of 
Eq. (44). The second-order term proportional to (T e + Ti) 
originates from perturbations of the particle fluid pressure. 
To the leading-order accuracy, the term 

(p-GQ(p-Ge)(Te+TQ m (fc-GaXfe-GelXTe+T,) 



p • (G e — Gi)e n, ■ \u e l — 

_ ikj [! + (! + «g)fcf Al] [! + (! + ^)k\/kl] (T e + Tj 
(/t e + Ki)(l + i>k)eB 

(50) 

contributes directly to jg. The combination of the last 
two terms describes the major linear wave dissipation due 
to particle diffusion and recombination. For fully magne- 
tized electrons, K e = 0. e jv a 1, and unmagnetized ions, 
Ki = fli/fi <§C 1, Eq. (50) reduces to the conventional loss 
term, 



(Gei — G,i)e 

_2\ 7.2 , 



(p-Gi)(p-Ge)(r e + rQ 

p • (Ge - Gi)e 



ij>%klC 2 s 



where in this limit i/>g « (1 + /tefcf /fcx^e^i/^e^i)- 

The first term in the RHS of Eq. (44) includes the major 
instability drivers. Unlike the two last terms, to the leading 
accuracy it is a first-order term. To retrieve second-order 
corrections, we have to linearize it with respect to small 
perturbations oc i&r^ / v e,i and iV. Denoting the corre- 
sponding linear corrections by S(- • • ), we have 



(p-G e )(K ■ V i0 )-(p-Gi)(K -V e0 ) 



(k ■ V i0 )S 



P- (G e - Gi 
P-Ge 



- (k ■ V e0 )S 



P-Gi 



+ 



P- (G e - G 

(fc ■ G el )(sk ■ v l0 ) - (k ■ G el )(SK ■ Vip) 

k ■ (Gei — Gii) 



P ■ (G e - Gi) 



(51) 



According to Eq. (35), in the last term of Eq. (51) we have 
8k — — iV. Presuming a uniform magnetic field, B — Bb, 
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expressing V e ,io in terms of Uo according to Eqs. (32), (33), 
and using Eqs. (47), (48), we obtain 



(fc • G el )(8K ■ V i0 ) - (k ■ G el )(6K ■ V i0 ) 



k-(G el -G n ) 

(/tj — K e )V • Uo 
1 +1pk 



+ K e Kib ■ (V X U ) 



(52) 



Recall that V • Uo can be expressed in terms of Vno accord- 
ing to Eq. (31). In this sense, the term proportional to V-Uo 
can be considered as an additional small contributor to the 
CD instability driving. The last vortex term oc V x Uo is un- 
related to any density gradients, but it may also contribute 
to 7s . 

Next we turn to the combination of the two first terms in 
the RHS of Eq. (51). The difference between the two frac- 
tions p ■ G e J p ■ (G e — Gi) and p ■ Gi j p ■ (G e — Gi) equals 
a constant value of 1, so that their linear perturbations are 
equal and we have 



(k ■ Vw) S 

- (k-Veo) 8 
where 

13 = 8 



P-G e 



P- (G e - Gi) 
P-Gi 



P- (G e - Gi) 



= —P(k ■ Uo), 



P-Ge 


= 8 


P-G t 


p- (Ge-Gi)_ 




/•(Ge-Gi). 



(k ■ Gel)(fc ■ SGj) - (k ■ Gil)(k ■ SG e ) 

[k ■ (Gei - Ga)] 2 
(k ■ G e i) (Ga -8p)-(k- G tl ){G el ■ 8p) 
[k ■ (Gei - Gn)] 2 



(1) 

Pfb + /3gd, 
(54a) 
(54b) 

(54c) 



Here the coefficient /?fb includes particle-inertia contribu- 
tions responsible for the FB instability, while /3gd includes 
the background density gradients responsible for the CD in- 
stability. 

First, we calculate the FB-instability driving term, 
— /? FB (fc ■ Uo). According to Eqs. (40) and Eq. (41), we 
obtain 



k-8G e 



k-8Gi 



B 



1 - K 



+ 



Hj% k _l 

~B~ 



(i + k 2 ) 2 ' k\ 

2 fc 2 
k 2 



n (e) 



1-Kt 



2\2 



fci 



Using Eq. (49), we have k ■ G e ,n, 

k- (Gei- Ga) n£ 



k ■ G e 



e) 
fci 



k-Uo 



k- (Ge 



k ■ Gn = - 



where according to Eq. (48), 



Gn) 0£ 



(0 



k-Uo 



r , =J =; . iKiK e ( Ke + Ki) (1 + 1p % ) k\ 

k ' (Gel " Gll) - + + 



As a result, we obtain from Eq. (54b) 
-j3 FB {k-U ) 



(H e + Ki) (1 + 
(1- K 2 )(1+ K 2 ) 



(1 + K 2 )Ke.Vi 



| (1- K 2 )(1 + K 2 ) | 
(1 + K 2 )KiV e 



(55) 



fc 2 



Here the first and second terms in the square bracket orig- 
inate from the ion and electron inertia, respectively. These 
terms correspond to particle oscillation in the perpendicu- 
lar to Bo plane. Notice that for K e > 1, the electron inertia 
terms becomes negative, meaning that it opposes the FB 
instability. The physical nature of this electron-inertia sta- 
bilization is fully analogous to that for the ion inertia above 
the magnetization boundary, Ki > 1 [Dimant and Oppen- 
heim, 2004]. The last term in the square bracket of Eq. (55) 
includes the effects of the electron {n e /i> e ) and ion (ni/vi) 
inertia in the parallel to Bo direction. The second and third 
terms in the square bracket are usually neglected, but the 
conditions for such neglect have not being properly analyzed 
and justified. We do this below. 

Similarly, we calculate the term — /3gd(& • Uo) that de- 
scribes the local gradient drift instability or stabilization. 
According to Eq. (37), 8p — —i Vno /no, and we obtain 
from Eq. (54c), 



-/3 G n(k-U ) 



'< < («e — Ki) ( fc|| • — ° — k± 



i(l + K 2 e )(l + K 2 )(k ■ Up) 
K e Ki(K e + Ki){l + l/>fc) 2 fc 2 

V ± n k\ 



no 



l + Cl + KiKe)^ 



(k X 6) • 



no 

V±n 
n 



fc 2 



(56) 



To summarize, we obtain for the combined FB and CD 
linear growth rate the general expression 



where 



7fb = — 



= 7fb + 7gd — jt — 2ano, 
1 + k 2 c ) 



(57) 



fci fci 



K e (1 + Vfe) 

(1-k 2 ) (1 + k: 

(1 + K 2 ) KiV e 



(1 + K 2 ) K^Vi 



+ 



K e ^ Ki 
V e Vi 



k 2 
k 2 



(58) 



7GD 



{1 + K 2 e ){l + K 2 ){k ■ Up) 
K e Ki(K e + Ki)(l + 1p k ) 2 k\ 



Vyno r V±n k 

X < (K e — Ki) fell • k± ■ —. 

* ' ' n n k- 

V±n 



l + (l + KiK e ) ^ 



(k x 6) • 



no 



(59) 



fcj [! + (! + Kpkf/kj] [! + (! + ^ 2 )fc 2 /fcj] (T e + Tj) 

(K e + Ki)(l + ^)eB 



7 T = 

\K e -r ni)\L -r wt)t-D 

(60) 

and the first-order shifted frequencies are given by 

Eqs. (36) and (39) from Dimant and Oppenheim [2011]. 

For magnetized electrons at a sufficiently high altitude 
where K e 1 and tp± = (fCe^) -1 <ti Ve/fi — 10 hold to- 
gether (see below), presuming |fc|||/fcx < k^ 1 k^ 1 , wc 
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obtain 



fei 



k-Up 
1 + 



so that 



7fb 



fci 



^(i-^)(^ ) 2 



7T 



(61) 



(62) 



7GD 



(i + k; 



)n ( , 4) 



Ki(l + ^jE)*i 



Vyno 
n 



(kxb) 



n 



(63) 

where = (T e +Ti)/rrii. Equation (62) shows that for 
fluid particles ion inertia is a destabilizing factor (7fb > 0) 
for all altitudes below the magnetization boundary defined 
by Ki = 1. For higher altitudes with m > 1, ion inertia is a 
stabilizing factor, 7fb < 0. The physical nature of this pe- 
culiar feature has been discussed in Dimant and Oppenheim 
[2004]. Note that the multipliers (l + ip%) in Eqs. (61) to 
(63) differ from the corresponding multipliers in the well- 
known expressions [Fejer et al, 1984] with the traditional 
definition of ip by a term tp±Ki — Ki/K e = 8q ~ 1.8 x 10~ 4 , 
that can be neglected. 

Now we discuss some of the conventionally neglected 
terms in the linear dispersion relation. The second and 
third terms in the square bracket of Eq. (55) have never 
been taken into account. In the case under consideration, 
the ratio of the second term to the first one is Vi/ '(mKeVe) = 
ip±(vi/v e ). In the lower ionosphere, we have z/j/i'e — 0.1, so 
that the electron inertia in the perpendicular to Bq plane 
can be neglected provided ip± <g; 10. This condition is ful- 
filled well above the 90 km altitude, i.e., in essentially within 
the entire altitude range of a possible FB instability devel- 
opment. At lower altitudes (the upper D region), the elec- 
tron inertia can become an additional stabilizing factor for 
the electron thermal driven instability [Dimant and Sudan, 
1997]. 

For the third term in the square brackets of Eq. (55), 
the ratio of its two sub-terms is (K e /ve)/{Ki/vi) = 
(mi/m e )(ui/u e ) 2 ~ 550, so that in the parallel to Bo di- 
rection it is the electron inertia term oc K e /v e that largely 
dominates. Then the ratio of the third term to the first one 
becomes (v e / Vi)^ / k\ ~ 10(fc||/fc±). It is of order unity or 
larger for the wave aspect angles, 9 = arctan(fc||/fcjj > 20°. 
Thus, particle inertia in the parallel to B direction, which 
is largely due to electrons, can only affect waves whose 
wavevectors lie well off the perpendicular plane, usually out- 
side the angular domain of linear instabilities. Such waves 
are heavily damped and hardly play a role in any linear or 
non-linear processes. This means that for waves in the lower 
ionosphere the last term in the square bracket of Eq. (55) 
is probably never of importance. 
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